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ABSTRACT 

Using the Southeastern Association for Research in Astronomy 0.6 meter telescope located at Cerro 
Tololo, we searched for variable stars in the southern globular cluster NGC 6584. We obtained images 
during 8 nights between 28 May and 6 July of 2011. After processing the images, we used the image 
subtraction package ISIS developed by Alard (2000) to search for the variable stars. We identified 
a total of 69 variable stars in our 10x10 arcmin^ field, including 43 variables cataloged by Millis & 
Liller (1980) and 26 hereto unknown variables. In total, we classified 46 of the variables as type RRab, 
with a mean period of 0.56776 days, 15 as type RRc with a mean period of 0.30886 days, perhaps 
one lower amplitude type RRe, with a period of 0.26482 days, 4 eclipsing binaries, and 3 long period 
(P > 2 days) variable stars. As many as 15 of the RRab Lyrae stars exhibited the Blazhko Effect. 
Furthermore, the mean periods of the RR Lyrae types, the exhibited period/amplitude relationship, 
and the ratio of N,,/ {Nab + A^c) of 0.25 are consistent with an Oosterhoff Type I cluster. Here we 
present refined periods, V-band light curves, and classifications for each of the 69 variables, as well as 
a color-magnitude diagram of the cluster. 

Subject headings: stars: variables: general-Galaxy: globular clusters: individual: NGC 6584 



1. INTRODUCTION 

Globular clusters provide an ideal controlled environ- 
ment to study variable stars and hence stellar evolution 
since the stars in the cluster presumably formed out of 
the same cloud of gas and dust, have the same compo- 
sition, have the same distance, and have the same age 
(Smith 2004). In addition to understanding the impor- 
tant aspects of late-stage stellar evolution in a controlled 
environment, variable stars in globular clusters are one 
of the primary standard candles that are used to measure 
the distances to globular clusters in the Milky Way and 
nearby galaxies. They also provide a method to easily 
obsca've thousands of stars at once, a useful characteris- 
tic when searching for variable stars. 

NGC 6584 is a relatively understudied globular cluster 
located 13.5 Kpc from the Sun in the southern constel- 
lation Telescopium (Harris 1996). It is an Oosterhoff 
type I cluster (Oosterhoff 1939) and Sarajedini & For- 
rester (1995) have estimated its metallicity at [Fe/H]=- 
1.44±0.16. Early studies of this cluster were performed 
by Bailey (1924) who identified nine potential variable 
stars in its vicinity. Millis & Liller (1980) published the 
locations of 48 variable stars in NGC 6584. They uti- 
lized photographic plates in their observations and anal- 
yses. As a result, the light curves that they present are 
relatively incomplete and the periods that they provide 
are somewhat imprecise. Sarajedini & Forrester (1995) 
confirmed these 42 RR Lyrae variables and tentatively 
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identified 56 possible yet undiscovered RR Lyraes by an- 
alyzing the variation in the frame-to-frame dispersion in 
V and B-V as a function of V for the stars in the in- 
stability strip of their color-magnitude diagram (CMD). 
Those stars which showed a large range in V dispersion 
were identified as possible RR Lyraes. In the same year, 
Samus et al. (1995), also using photographic photome- 
try, again confirmed Millis & Liller's 48 variables. Since 
the results of Sarajedini & Forrester (1995) and Samus et 
al. (1995) image subtraction software and better CCDs 
have become available. 

With the combination of CCD images, an observation 
window spanning more than 5 weeks, and an image sub- 
traction method developed by Alard (2000), we are able 
to detect and analyze hereto undiscovered variable stars 
and produce unprecedented high-quality light curves and 
periods for each variable. In this study we use Alard's 
image subtraction method to search the central 10x10 
arcmin^ of globular cluster NGC 6485 for variable stars 
from observations obtained in May, June, and July of 
2011. We present a preliminary color-magnitude diagram 
for the cluster, as well as complete light curves and peri- 
ods for 43 of Millis & Liller's variables and for 26 newly 
identified variables. 

2. OBSERVATIONS AND REDUCTION 

We obtained images using the Southeastern Associ- 
ation for Research in Astronomy (SARA) 0.6 meter 
telescope at Cerro Tololo Interamerican Observatory 
(CTIO). Observations were made on the nights of 28 
May; 2, 22, 25, and 28 June; 2, 5, and 6 July 2011 us- 
ing an Apogee Alta E6 camera with a 1024x1024 pixel 
Kodak KAFIOOIE chip, a gain of 1.5 electron per ADU 
count, and an RMS noise of 8.9 electrons. We used 1x1 
binning resulting in an image scale of 0.61"/px resulting 
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Fig. 1. — The inner 5x5 arcmin^ of the ref .fits image used for 
all nights of the analysis/image subtraction. This image is a com- 
bination of the best seeing images, with the combined seeing near 
1.3". Note that the 6 brightest stars in the image are overexposed. 

in a 10x10 arcmin^ field of view. The CCD was kept 
at a constant temperature of -25 Celsius. The typical 
seeing was 1.7" and ranged from 1.2 to 2.8". On all 
nights of observation, 240 second exposures were taken 
using a Bcsscl V filter. The exposure time allowed us to 
detect dimmer variables but caused 6 of the stars in the 
field many to be overexposed. Most of these overexposed 
stars appeared to be in the foreground of the cluster. On 
the nights of 2, 5, and 6 July 2011, additional exposures 
were taken using Bcsscl B, R, and I filters with exposure 
times of 300, 90, and 50 seconds respectively. All images 
were debiascd, flat-fieldeci, dark subtracted, and cleaned 
of hot, cold, and bad pixels using MaximDL. 

3. ANALYSIS 

Image Subtraction was accomplished using the ISIS- 
2.2 package (Allard & Lupton 1998; Allard 2000). ISIS 
accounts for changes in seeing conditions by convolving 
a high quality reference image to the point spread func- 
tions of each individual frame and then subtracting the 
two images to determine any change in intensity. The 
package contains seven c-shell scripts and three parame- 
ter files that were used in our analysis. The different c- 
shell scripts are dates, interp (registration and interpo- 
lation) , ref (build reference frame) , subtract (subtracts 
images from reference), detect (stacks subtracted im- 
ages), find (finds variables above a user-defined thresh- 
old), and phot (produces light curves for all detected 
possible variables), along with three parameter files. 

We began by running a modified version of the dates 
script, which identified and registered each image file 
in the image directory. We then ran the interp script 
which aligned each image to a previously chosen align- 
ment image. Next, the ref script combined five of our 
highest quality images into a reference frame (Fig. 1). 
In order to maintain consistent and comparable values 
for relative flux across nights, one reference image was 



Fig. 2. — The inner 5x5 arcmin^ of the var.fits frame ouput 
from ISIS for all eight nights from SARA CTIO. The relative 
amount of variation is indicated by the brightness of the star. The 
6 overexposed stars show significant variation only because ISIS is 
unable to fit an accurate point spread function to them. 

created using five of the best images (those with see- 
ing between 1.2 and 1.3"). We then ran the subtract 
script, which first convolves the high-quality reference 
frame to the point-spread function of the image to be sub- 
tracted (in order to compensate for variations in seeing 
over the course of the observations), and then subtracts 
the convolved reference frame from each image. These 
subtracted images were then stacked using the detect 
script creating a var.fits image which shows the de- 
gree of variability across the observation run for each 
star in the cluster as shown in Fig. 2. We then used 
the find script to identify variable stars above a cer- 
tain user-defined threshold called SIGTHRESH. We used a 
SIGTHRESH of 0.03 in our analysis to ensure that all pos- 
sible variables were found. This value typically resulted 
in 1500 possible variables. Finally, the phot script was 
used to generate light curve data files for each identified 
variable. These light curves allowed us to eliminate false 
positives from the find script. Variable stars produce 
very distinct light curves, and thus false positives due 
to noise and variations in seeing were easily identified 
amongst these plots and discarded. 

Periods were determined, when possible, based on com- 
bined photometric data from all eight nights. Wc used 
the period finding software AVE (Analisis de Variabili- 
dad Estelar, Analysis of Estcllar Variability) from Grup 
d'Estudis Astronomies. AVE uses the Lomb-Scargle al- 
gorithm (Scargle 1982) to determine the most probable 
period for a given set of data. Due to the more than 5 
week observation span, wc were able to determine the pe- 
riods of most of the variables with periods less than 1 day 
to an accuracy of 10~^ days. Those variables exhibiting 
the Blazhko Efi'ect were less accurate due to modulation 
of their light curves. 

Millis & Liller's (1980) variables were identified visu- 
ally and astrometrically and then confirmed by compar- 
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TABLE 1 

Classification Summary 



Variable Type 


Count 


Period (days) 


RRab 


46 


0.56776 


RRc 


15 


0.30886 


RRe 


1 


0.26482 


Eclipsing 


4 


0.32431 


LPV 


3 


> 20 



ing their periods with our own. Five of the variables 
identified by Millis & Liller were outside the field of view 
of our images thus were not be observed. 

4. RESULTS 

Table 1 summarizes our findings by classification. Ta- 
ble 2 list the positions, periods, and classification of all 
of the variables that we identified. In Table 2 variables 
previously identified by Millis & Liller (1980) are labled 
with a V and the number they assigned in the first col- 
umn. Newly discovered variables are labeled as NV and 
are in order of increasing right ascension. We detected 
all the Millis & Liller (1980) variables that were within 
our 10x10 arcmin^ field. Their variables 24, 27, 40, 43, 
and 44 were outside our field of view. We did note that 
variable 45 was misidentified on the finder chart in Mil- 
lis & Liller (1980), thus the position we list for of this 
variable has been shifted by several arc-seconds in Ta- 
ble 2. We found periods for variables V12, V19, and 39 
that Millis & Liller were unable to determine. We also 
were able to eliminate those periods where multiple pos- 
sible periods were listed; specifically V5, V20 and V35. 
Besides the increased accuracy of our periods the only 
significant change from Millis and Liller was V13 which 
we found to be 0.34549 days as compared to 0.39 days. 

Each detected variable was classified based on the 
shape of its light curve, its period, and its maximum 
amplitude. In addition to the Miller & Liller variables 
we found 26 additional variables in our field, most of 
which are RR Lyrae variables, long period variables, or 
eclipsing variables. Of the total 69 variables detected, 46 
were classified as RRab stars, 15 as RRc, 1 possible RRe, 
4 eclipsing binaries, and 3 long period variables. From 
the fraction Nc/{Nab + Nc) we confirm that the cluster 
is of Oosterhoff Type I. RRab stars are characterized by 
large amplitudes and periods between 0.4 and 0.8 days. 
RRc stars have smaller amplitudes and periods between 
0.2 and 0.4. NV9 has the characteristics of a type RRe 
variable (Alcock et al. 1996; Bono et al. 1997). It has a 
lower amplitude and period than type RRc stars. Eclips- 
ing binary stars are easily identified by their unique light 
curves. As many as 15 of the RRab stars appeared to ex- 
hibit the so-called Blazhko effect (Blazhko 1907; Smith 
2004), which is characterized by a long term modula- 
tion in amplitude and period. These variables are noted 
with and asterisk in Table 2. Their periods are likely 
to be less accurate given their behaviour. We have not 
assigned periods to the long period variables at the tip 
of the asymptotic giant branch because further observa- 
tions will be needed to establish the their periods. We 
assign 'LPV to these variables with a period greater than 
2 days. 



TABLE 2 
Variables Stars in NGC 6584 



v# 


RA (h,: 


m,s) 


Dec (°,' 


, ) 


Period (d) 


Type 


VI 


18 


18 


28.82 


-52 


13 


27.0 


0.49520 


RRab 


V2 


18 


18 


42.20 


-52 


12 


30.0 


n 54465 


RRab 


^V3 


18 


18 


42.14 


-52 


13 


16.2 


0.47162 


RRab 


V4 


18 


18 


40.30 


-52 


13 


25.5 


0.54377 


RRab 


V o 


18 


18 


40 49 


-52 


13 


9Q 1 
zy . _L 


n 9079^^ 
u.zy i zo 


RRc 


V u 


18 


18 




-52 


12 


7 R 


u . uzu 1 y 


r\,r\,ciij 


V 1 


18 


19 




-52 


11 


OO.O 


^A7A^ 

U.04 1 '-to 


IXlXcilj 


V o 


18 


18 




-52 


13 


OO.O 


U.OUO-L -L 


RRc 


ayg 


18 


18 


43.55 


-52 


12 


2.7 


U.4:yO -LO 


RRab 


VIO 


18 


18 


37.07 


-52 


11 


36.7 


0.69432 


RRab 


ayil 


18 


18 


35.42 


-52 


13 


10.8 


0.53945 


RRab 


»V12 


18 


18 


35.84 


-52 


12 


32.6 


0.59505 


RRab 


V13 


18 


18 


35.85 


-52 


14 


32.3 


0.34549 


RRc 


avi4 


18 


18 


41.20 


-52 


13 


43.2 


0.59954 


RRab 


V15 


18 


18 


23.56 


-52 


12 


41.9 


0.63632 


RRab 


V16 


18 


18 


42.04 


-52 


12 


58.6 


0.58220 


RRab 


V17 


18 


12 


41.87 


-52 


13 


22.3 


0.56338 


RRab 


^V18 


18 


18 


32.99 


-52 


13 


18.3 


0.55789 
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V23 


18 


18 


30.76 


-52 


13 


10.0 


0.59812 


RRab 


V25 


18 


18 


56.91 


-52 


11 


12.8 


45336 


RRab 


V26 


18 


18 


33.62 


-52 


13 


52.1 


0.56970 


RRab 


V28 


18 


18 


48.09 


-52 


12 


19.8 


0.49815 


RRab 


V29 


18 


18 


41.49 


-52 


11 


40.5 
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RRab 
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18 
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42.85 
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31.9 
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RRab 


V31 


18 
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-52 


13 


38.5 
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RRab 


V32 


18 


18 


40.71 


-52 


14 


5.4 
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RRab 


V33 


18 


18 


26.81 


-52 
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0.28595 


RRc 
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13 
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^'NVig 


18 
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-52 
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NV20 
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8 
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EB 
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EB 



variable cxhibitiing Blazklio Effect. 
^ variables with significantly revised periods, 
misidentified by Millis & Liller (1980). 




Fig. 3. — Phased light curves for identifieded variables with periods under a day. The vast majority of these variables are RR Lyrae 
stars. 
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Fig. 3. — (Continued) Phased light curves for identified variables with periods under a day. The vast majority of these variables are RR 
Lyrae stars. 
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Fig. 3. — (Continued) Phased light curves for identified variables with periods under a day. The vast majority of these variables are RR 
Lyrae stars. 
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Fig. 3. — (Continued) Phased light curves for identified variables with periods under a day. The vast majority of these variables are RR 
Lyrae stars. 
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Fig. 3. — (Continued) Phased light curves for identified variables with periods under a day. The vast majority of these variables are RR 
Lyrae stars. 
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Fig. 4. — Light curves for the 3 long period variables found in our analysis. 



Cluster membership was determined using a color mag- 
nitude diagram (CMD). The CMD was constructed from 
137 V and 155 R frames from several nights to determine 
magnitudes, colors, and cluster membership of the vari- 
ables. DAOPHOT was used to determine V and R magni- 
tudes of the stars (Stetson 1987). In the CMD shown 
in Figure 5 RRab and RRc variables are indicate by the 
red and blue colored points, respectively. RRc variables 
tend to lie on the warmer (left) side of the instability 
strip whereas RRab variables are to the cooler (right) of 
the instability strip. Scatter in the variables is partly 
due to blending of the variables with other stars, our av- 
eraging of several hundred images before using DAOPHOT, 
and possible variable reddening in the direction of the 
cluster. 

Six of our variables, Vll, V22, NV7, NV19, NV20, and 
NV21 were blended with other stars and show significant 
offset to the right of the instability strip and are likely 
cluster members. Four of the RRab variables V15, V41, 
NV4, and NV15 were significantly brighter than the RR 
Lyrae variables found in the instability strip of the cluster 
and are suspected not to be cluster members. From the 



CMD we are able to determine that 2 of the eclipsing 
variables to the right of the giant branch are not cluster 
members. The 3 long period variables that we found 
all lie near the tip of the asymptotic giant branch. The 
remainder of the RR Lyrae stars lie very near or in the 
instability strip and thus are likely cluster members. 

5. CONCLUSIONS 

Over a 39 day time span from late May to early July 
of 2011 we observed the globular cluster NGC 6584 for 
8 nights using the SARA 0.6 meter telescope at CTIO. 
Using the image subtraction software ISIS we were able 
to find precise periods of previously known variables and 
26 newly discovered variables stars, most of which are RR 
Lyrae stars. Within our 10x10 arcmin^ field of view we 
detected all the variables fisted by Milfis & Lillcr (1980) 
bringing the total to 74 variables in the vicinity of the 
cluster. Of all the RR Lyrae variables all but 4 appeared 
to be members of the cluster. Two of the 4 eclipsing 
binaries are likely to be cluster members. In total, we 
classified 46 of the variables as type RRab, with a mean 
period of 0.56776 days, 15 as type RRc with a mean 
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Fig. 5. — Color-magnitude diagram of NGC 6584, with variable stars indicated. Note that the RRc (blue) variables are bluer than the 
RRab's (red). All but 6 of the variables found appear to be cluster members. The two eclipsing variables to the right of the giant branch 
of the cluster are likely not cluster members. The 2 RRc and 4 RRab Lyrae variables just to the right of the instability strip are all blends. 
The 4 brightest RRab's may be foreground objects. 



period of 0.30886 days, perhaps one lower amplitude type 
RRe, with a period of 0.26482 days, 4 eclipsing binaries, 
with a mean period of 0.3243 days, and 3 long period 
(P > 2 days) variable stars. As many as 15 of the RRab 
Lyrae stars exhibited the Blazhko Effect. The ratio of 
Nc/{Nab + Nc) of 0.25 is consistent with cluster being of 
type Oosterhoff Type I. 
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